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We examine the dephasing rate of a Josephson charge qubit system due to background 
charge fluctuations. We consider single qubit and two-charge traps. The transition probability 
was controlled to a state where two traps were occupied. The transition probability was 
affected by the Coulomb blockade effect that occurs between two charge traps. To obtain 
the dephasing rate, we computed the spectra of random frequency modulation signals. Our 
results show that the interaction between charge traps suppresses dephasing. 
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Among various proposals, quantum bits (qubits) in solid state materials, such as, super- 
conducting Josephson junctions 1 ) and quantum dots, 2-4 ) have the advantage of scalability. 
Proposals to implement a quantum computer using superconducting nanocircuits are prov- 
ing to be very promising 5 ' 6 ) and several experiments have already highlighted the quantum 
properties of these devices. 7 ) Such a coherent- two-level system constitutes a qubit and the 
quantum computation can be carried out as the unitary operation functioning on the multiple 
qubit system. Essentially, this quantum coherence must be maintained during computation. 
However, dephasing is hard to avoid due to the system's interaction with the environment. In 
terms of a bonding-antibonding bases, the decay of off-diagonal elements of the qubit density 
matrix signals that dephasing is occurring. This dephasing is characterized by the dephasing 
time T2. Various environments can cause dephasing. In superconducting nanocircuits various 
sources of decoherence are present, 5 ) such as fluctuations originating from the surrounding 
circuit, quasiparticle tunneling, background charge fluctuation (BCF), and flux noise. For a 
charge qubit system, BCF is one of the most critical dephasing channels. 8 ~ n ) 

BCFs have been observed in various kinds of systems. 12 ~ 15 ) In nanoscale systems, they are 
the electrostatic potential fluctuations due to the dynamics of electrons, or holes on a charge 
trap. In particular, the charge at a charge trap fluctuates with the Lorentzian spectrum form, 
which is called random telegraph noise in the time domain. 14,16 ) The random distribution of 
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the positions of such dynamical charge traps and their time constants lead to BCFs or 1/f 
noise. 17 ) In solid-state charge qubits, these BCFs result in a dynamical electrostatic distur- 
bance and hence, dephasing. The theoretical effect of 1/f noise on a charge Josephson qubit 
has been examined previously. 8 ~ n ) 

We investigated how the electrostatic disturbance coming from two or more dynamical 
charge traps affects the quantum coherence of a qubit. In past studies, an environment com- 
posed of free charge traps had been considered. 8 ' 10 ) When such an environment is interacting 
with itself, its characteristic nature would be expected to affect the relaxation phenomena. 
In present study, we especially concentrated on the correlation effect between the charges in 
the traps. We consider pure dephasing as an event which occurs when the dynamical charge 
traps induce fluctuation in extra bistable bias. It should be noted that this dephasing process 
does not mean the qubit is entangled with the environment, but rather, that the stochastical 
evolution of an external classical field, is suppressing the off-diagonal density matrix elements 
of the qubit after being averaged out over statistically distributed samples. 

The system under consideration is Cooper pair box. 5 ) Under appropriate conditions ( 
charging energy Eq much larger than the Josephson coupling Ej and temperatures ksT -C 
Ej) only two charge states are important, and the Hamiltonian of the qubit H q b reads 

u 5Ec ± Ej m 

where the charge bases {|0 >,|1 >} is expressed using the Pauli matrices, and the bias 
5Eq = Ec{l — C x V x /e) can be turned by varying the applied gate voltage V x . The environment 
is a set of BCF electrostatistically coupled to the qubit, 8 ' 10 ' 18 ' 19 ) 

H qb _ imp = ^2^^a z (4di-^) (2) 
i=i 

where d\ and di are the electron creation and annihilation operators of a charge trap, i is the 
index of N charge traps, and the coupling with the qubit is such that each BCF produces a 
bistable extra bias hJci- Because qubit Hamiltonian consists of Ej and 8Eq, the dephasing 
consists of that with dissipation and pure dephasing. In general, the dephasing with dissipation 
can be neglected as follows. For physical setups, 5Ec — 122 /ueV, and Ej ~ 34 //eV; 9 ) By 
perturbation method, 5 ) the ratio of the dephasing rate with dissipation to pure dephasing rate 

E 2 \2 

is roughly given by ^gr ^ E -i + ^ )/n 2 +X 2 111 ^ ne P resence °f the bistable extra bias, where A is 
the transition rate of the dynamical charge trap. For the above experimental setups with the 
dynamical charge trap with low frequency, we can neglect the effect of Ej because Ej < 5Eq 
and v c h J > A. Then the pure dephasing event is critical. In final results of present study, 
we discuss about the many charge traps which are interacting with each other. For this case, 
the dominant process is different, we discuss about this behavior latter. We neglect the back 
action from the qubit to charge traps, namely, the transition rates of charge traps do not 
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depend on the qubit state. This assumption is justified because the qubit induces static shift 
of the chemical potential of charge traps. And the change in the chemical potential of a charge 
trap does not renormalize the transition rate of itself. 20 ) 

Using the environment variable Xj(t)(=< d\{t)di{t) > r —1/2), where < A(t) > r is trace of 
the operator A(t) about the electron reservoir of the charge trap, we rewrite the Hamiltonian 
in terms of the Pauli matrix as 

H= 6 ^a z + f:^a z X l (t), (3) 

i 

we assume that the charge traps are strongly coupled with their charge reservoirs and the 
time evolution of Aj(t) is a Poisson process. 

Following the time evolution of density matrix of qubit, we obtain the following off- 
diagonal element pu{t) = p\2{to)e l&Ec ^ h ^ t ~ t °^ e ■^ t » dTX ^ T \ where t$ is the initial time, and 
x(t) = 'Yld JaXi{t) takes 2^ with possible different values of ai,...,a 2 iv. The fluctuation 
in tunneling coupling constant is pure dephasing and does not accompany relaxation of the 
population. Therefore, diagonal elements of the qubit density matrix do not change. 

In the following, we estimate the ensemble average of off diagonal element of the density 
matrix, E[p\2{t)\ = pi2(to)e lSEc VK*-*") < e l -h dTX ( T ) > _ p or this q Uan tity, we can apply 
the characteristic functional method, 21 ' 22 ) namely, R(t) =< e~ l $a drx ( T ) >= ^ 2 fc=1 Pi-Rjfc(i), 
where pi is the occupation probability of the state i, which can be determined by the stationary 
condition, = —piPi{t) + X^=i ^i^jiPj{t)> w here \j is the transition probability defined for 
very short time At with i ^ j, and pi = Ylj^i ^ij ^ s the emission rate during this time. The pi 
has the properties YliPi = 1 an d PiPi = ^jiPj- The average of x, r], is given by YliPi a i 



and the variance a is given by yYliPi a i ~ V 2 - The function Rik satisfies following the first 
order differential equation, 

- = {idk ~ Pk)Rik(t) + ^2 X mk R km (t) 

m^k 

= ^RimAmk, (4) 

m 

with the initial condition Rij(t = 0) = 5{j. T^ 1 characterizes the exponential tail of long-time 
dephasing behavior. This quantity is obtained by Min(— Re(ej)), where ej's are the eigenvalues 
of A. While for very short t, the R(t) shows Gaussian behavior. 8 ) 

First, we examine the single charge trap case (N=l) where we have 

A= fi ai -x u X u \ (5) 

V ia.2 - \d J 

where the X u (Xd) is the transition rate from the 1st state to the 2nd state ( 2nd state 
to the 1st state). The resultant T 2 -1 is given by T^ 1 = \{X U + Xd — ReVA) where A = 
4aid2 + 4i(aiArf + d2X u ) + (X u + Xd — ia\ — id2) 2 , and Jq = a.2 — «i- For weak coupling 
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Fig. 1. Scheme of two charge trap system. The wells 
represent the charge traps and the arrows represent 
the transition between each of four states. Ay's in- 
dicate the transition probabilities from state i to 
state j. 



(\Jc\ ^ max(A u , Xd)), Jc = d = —2a\ = 2a>2, and d characterize strength of bistable extra 
bias, T 2 is given by 

T 2 X = X u X d d 2 /(X u + X d f. (6) 

For strong coupling, (\d\ 3> X Ujd ), T 2 _1 = (A n + X d )/2. These results coincide with those found 
by Itakura and Tokura, 8 ^ where the dephasing time was derived using a different method. 

Next, we examined the two traps, (N=2) including the Coulomb blockade effect occurring 
between the traps. When two traps are located close to one other, there should be capacitance 
coupling between two occupied traps. However, we neglect tunneling between the charge traps. 
There are four states: both charge traps empty, left charge trap occupied, right charge trap 
occupied, and both charge traps occupied (Fig. 1). X'^s are transition rates from i state to j 
state, and we neglect the transition processes between 1 state and 4 state, and 2 state and 3 
state. In general, we notice A12 > A34 and A13 > A24, where the equations are for the absence 
of Coulomb blockade effect. 

For actual calculation, we restricted the parameters for the transition rate which are 
symmetric for two charge traps, A n = A12 = Xi^,X d = A21 = A3i,A^ = A24 = A34,A^ = 

A A' A A' A A' 

A42 = A43. The occupation probabilities are p\ = — ^-,P2 = P3 = -77^,^4 = U D U , where 
D = X d X' d + (2X' d + X' U )X U . 

First we discuss the high temperature behavior, where energy differences of each state 
are lower than temperatures, while 5Eq and gap energy of cooper pair of qubit are much 
higher than temperatures. To demonstrate the effect of Coulomb interaction transparently, 
we chose parameters, X u = X d = X' d = X and A^ = A', and calculate T2 while changing A' in 
the range < A' < A. Here, the occupations are pi = P2 = P3 = 3A ^ A / and p4 = 3A A hA , , and 
the amplitudes are, a\ = —d, (12 = 03 = and 04 = d. The dephasing rates for weak coupling 
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(d <C A) are 




£_ 

4A> 
2d 2 
27A' 



for A' = A 



for A' 



= . 



(7) 



In Fig. 2, the results by numerically solving Eq. (4) are plotted. This plot shows that two 
limits (TrT 1 = |^ and ) are smoothly connected in the intermediate parameter region 
for d/X = 0.1 (weak coupling case). In the limit of no interaction (A = A'), the dephasing 
rate becomes twice of that for the single charge trap. For a strong interacting limit (A' = 0), 
the time evolution reduces to that of a single charge trap with asymmetric transition rate 
of Xu U9le = |Ai and \ s ™ ale = |A X with Ai = |A and the dephasing rate is smaller than 
that of a single charge trap. For d/X = 2, there is a rapid increases in the dephasing rate 
when A ~ A'. There are four eigenvalues for characteristic equation of Ri m . Therefore there 
are four characteristic dephasing rates for this case. This singularity appears because two of 
them become same, thus the transition from weak coupling to strong coupling occurs there. 
Such a singularity also appears for the dephasing due to single charge trap. 8 ) All plots show 
that the dephasing rate increases with A'/A, which indicates that the effect of interaction, or, 
the screening effect, suppresses the dephasing compared with that of non-interacting charge 
traps. In this analysis, the average, r], changes with the ratio A'/A. If we choose such that 
the average of rj is invariant, the results are the same. The reason is that for both cases, 
«2 — o,\ = 03 — a\ = d and 04 — 02 = 04 — 03 = d, independent of A'/A, and the difference 
between the former case and the latter case only leads to the modulation of Rabi oscillation 
frequency. 

We also examined the Gaussian behavior, which is the short-time regime for t < 
mzn(^, 3/(max(A, A'))). 8 ** For dephasing due to a single charge trap, the off-diagonal element 
of density matrix decay is represented as, R(t) ~ exp(— i(^|-) 2 ), where T 2 ~ 2 is given by 
For two non-interacting charge traps, T^ 2 = and for strong interaction (A' = 0), T 2 ~ 2 = 
This behavior shows that the dephasing is suppressed as interaction increases, even for the 
Gaussian behavior. It should be noted that the decay rate of Gaussian behavior depends on 
the total charge of charge traps, where we chose zero as the mean of amplitudes. In present 
examinations, the effect of interaction between charge traps is discussed, while the numerical 
estimation of dephasing rate due to non-interacting BCF had been done in refs. 8-10. Please 
note that T 2 ~ 2 depends only on distribution of d 2 and number of charge traps, although T 2 _1 
depends on distributions of d 2 /X and number of charge traps; Eq. (7). The former result 
coincides with that of ref. 10 when the initial state is in thermal equilibrium. 

At lower temperatures than the energy differences of each state, we have asymmetric 
transition rates. Therefore, we must consider the effect of temperature. In order to satisfy the 
detailed balance condition, temperature and electron correlation leads to the following forms 

A A + E charge 

of the transition rate: 23 - 1 A^ = X' d = A, X u = Xe k B T and X' u = Ae k B T . The definitions of 
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Fig. 2. The A'/A dependence of dephasing rate 
T% X/cP- The lines indicate the numerical results 
where the parameters are d/X= 0.1, 2 and 4. 



energy difference are: A = E 2 - E 1 = E 3 - E u A + E charge = E 4 - E 2 = E A - E 3 . The E charge 
is the capacitive energy between two charge traps. In this case, the probabilities of population 
obey classical Boltzmann distribution, where pi = — 4 e ^wtu r • At high temperatures (T » 
Ei/kB), the occupation probabilities become p\ = p 2 = P3 = p A = 1/4. In Fig. 3, we show the 
dephasing rate obtained numerically with the amplitudes set to a\ = — d, a 2 = 03 = 0, 04 = d. 
We chose the numerical parameter of d/X = 0.1 (weak coupling). Using analytical expressions 
of dephasing rate for single charge trap, Eq. (6), the equation of the normalized dephasing 
rate for weak coupling case (d <C A) is given by T^ ingle = x ^^-^/fjr^ • The behavior 
of the traps for N=2 requires detailed examination. For weak interaction {E c h arge <C A), the 
dephasing rate due to the two charge traps is twice that of the dephasing rate due to the single 
charge traps. From Eq. (6), the dephasing rate becomes suppressed as the asymmetry of the 
transition rates increases. At low temperatures (ksT <C A), the dephasing rate is suppressed 
exponentially, because the asymmetry of the transition rates increases with a decrease in 
temperature. At high temperatures, (ksT » A), the dephasing rate is again suppressed. 
The reason is that, the characteristic transition rate (X u + A^), increases as the temperature 
increases. When coupling between the qubit and charge traps is weak (d <C max(A u , A^)), 
the magnitude of the fluctuations in the trace of a state on the Bloch sphere decreases with 
increasing, A u + A^, 8 ) hence the dephasing rate decreases as well. 
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Fig. 3. The A/ksT dependences of dephasing rate 
T^X/d 2 with d/\ = 0.1. The lines indicate 
the numerical results where the parameters are 
E char gel 'A = 0, 0.1, 10 and oo. The dot-dashed line 
indicates the analytical curve for single charge trap. 



Next, we must consider the behavior of the traps when Coulomb interaction is strong 
(Echarge ^ A). Except for very high temperatures (ksT 3> Echarge + A), the dephasing rate 
due to two dynamical charge traps is, T^X/d 2 = 2e~ A / fcsT / '{l+2e~^ kBT f . Then, comparing 
with Eq. (6), two charge traps are equivalent to a single charge trap with asymmetric transition 
rate, A^ = A, A„ = 2e -A / fcsT A. At intermediate temperatures (E c harge » fc^T » A), the 
dephasing rate is smaller than that of a single charge trap. The reason for this behavior is 
that two traps behave as a single charge trap with a larger characteristic transition rate, A(l + 
2 e -A/fc s T^ com p are d w ith that of the single charge trap A(l + e~ A / fcsT ). At low temperatures, 
(A, E c harge <C ksT), the dephasing rate decreases exponentially as the temperature decreases 
irrespective of E charge . 

Finally, we examine N identical charge traps which are located close to one another. To 
simplify the discussion, we consider a system of N charge traps symmetrically coupled with 
a qubit. There are: one empty state (i = 0), N single occupied states (i = 1), (N-l)N/2 two 
occupied states (i = 2). ■ ■ one fully occupied state (i = N). For strong and long-range 
Coulomb interaction, the empty state and single occupied states are relevant. When there 
is weak coupling, we have T 2 _1 = T n N ^ e -A/k B Tys > where d and A are coupling constants 
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between the qubit and the background charges, and characteristic transition rate of charge 
traps, respectively, which are identical for all charge traps. When there is non-interaction, N 
charge traps behave independently, and T 2 l = ^ ^^-A/tgTp ; hence the interaction between 
charge traps suppresses the dephasing rate. At high temperatures, the analytical solution of 

T- 1 and T' 2 are ffiven bv T ^ 1(str ° ng " ltoactio ") _ 8 ^(strong interaction) _ ± 

i 2 ana i 2g are given oy 27 1 (free) - (at+TF' 27/ (free) ~ iv> wnere 

we chose zero as the mean of amplitudes. Therefore, the dephasing rate becomes suppressed 
more effectively as the number of charge traps increases. It should be noted that when charge 
traps are interacting strongly each other, the dephasing rate with dissipation in the large N 
limit is given by ^2 ^_ E i j^r • Then the dephasing rate with dissipation becomes also suppressed 
with increasing N. While the dephasing rate with dissipation becomes gradually dominant over 
the pure dephasing rate as increasing N, we do not argue this effect since both rate vanish 
with N. 

In conclusion, we examined the dephasing rate of a two-level system, coupled with a 
classical environment made of N charge traps. The environment changes its bistable extra 
bias, which results in pure dephasing. When the charge traps fluctuate independently, the 
total dephasing rate is the simple summation of the dephasing rate of each charge trap. If 
multiple charge traps are interacting with each other, the dephasing rate is slowed, when T 
is not much smaller than A/fee- At high temperatures, (T 3> A/ks), more than one charge 
traps with large Coulomb interaction results in a smaller dephasing rate than that of the single 
charge trap. It should be noted that the other channels of dephasing exist, such a dephasing 
rate should be added to present dephasing rate. And present estimation of dephasing rates 
corresponds to that of free induction decay, 9 ) not that during gate operation in such a case 
the charge degeneracy state (SEq = 0) should be manipulated. The numerical evaluation of 
dephasing rate for such a situation has been done in refs. 8, 9. 
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